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ABSTRACT 

Chandra has detected optically thin, thermal X-ray emission with a size of ~ 1 arcsec 
and luminosity ~ 10 33 erg s _1 from the direction of the Galactic supermassive black 
hole (SMBH), Sgr A*. We suggest that a significant or even dominant fraction of 
this signal may be produced by several thousand late-type main-sequence stars that 
possibly hide in the central ~ 0.1 pc region of the Galaxy. As a result of tidal spin- 
ups caused by close encounters with other stars and stellar remnants, these stars 
should be rapidly rotating and hence have hot coronae, emitting copious amounts 
of X-ray emission with temperatures kT<^ a few keV. The Chandra data thus place 
an interesting upper limit on the space density of (currently unobservable) low-mass 
main-sequence stars near Sgr A*. This bound is close to and consistent with current 
constraints on the central stellar cusp provided by infrared observations. If coronally 
active stars do provide a significant fraction of the X-ray luminosity of Sgr A*, it should 
be moderately variable on hourly and daily time scales due to giant flares occurring 
on different stars. Another consequence is that the quiescent X-ray luminosity and 
accretion rate of the SMBH are yet lower than believed before. 

Key words: Galaxy: centre - stars: coronae - X-ray: stars. 



1 INTRODUCTION 

Our Galactic nucleus with its 
(SMBH) of mass M BH ~ 4 x 10° M e 
intensive studies in radio, submillimeter 
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black hole 
a subject of 
X-ray 



is 

infrared 



and g amma-ray bands fsee lGenzel. Eisenhauer. fc Gillessenl 
I2010J for a recent review). The Chandra X-Ray Observa- 
tory, with its excellent spatial resolution ( Weisskop f et al.l 
2002), has revealed a rich variety of compact and diffuse X- 
ray sources in the central square degree of th e Milky Way 
llMuno et alj 12004 IWang. Dong, fc LandliuOrj : iMuno et all 
l200Sl . 12009). One of the most interesting findings is the de- 
tection of faint extended X-ray emission surr ounding the ra- 
dioso urce Sgr A* associated with the SMBH (Ba ganoff et al.l 
12001 hereafter B03). The X -ray source has a radius of only 
~ 1", or ~ 0.04 pc for a Galactic Centre (GC) distance 
of ~ 8 kpc. Its X-ray luminosity is few 10 33 erg s _1 , and 
the spectrum is consistent with thermal emission from op- 
tically thin kT ~ 1-4 keV plasma, absorbed by a column 
Ah ~ 10 23 cm -2 of neutral gas (B03). This column density 
is roughly consistent with the interstellar extinction towards 
the GC region. 

The currently favoured view is that the X-ray emission 
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from the central arcsecond around Sgr A* is truly diffuse 
and ass ociated with ho t gas flowing towards the central 
SMBH (iQuataertl 120021: B03: lYuan. Quataert. fc Naravanl 
120031 : ICuadra et al.lbo'orj ; IXu et al.ll2006l h According to this 
scenario, the X-ray luminosity is dominated by thermal gas 
emission near the Bondi radius, _Rb = GMbh/c 2 ~ 0.05 pc 
(where c s is the gas sound speed, corresponding to a tem- 
perature of a few 10 7 K). The Chandra observations im- 
ply that the outer inflow rate is M ou t ~ 10 _5 Mq yr _1 . 
Although this rate is very low compared to active galac- 
tic nuclei, if the SMBH were accreting gas at this rate in 
radiatively efficient mode its bolometric luminosity would 
be several xlO 40 erg s _1 , i.e. 5 orders of magnitude higher 
than the total luminosity of Sgr A* (predominantly emit- 
ted at radio to submillimeter wavelengths). The radiative 
output can be much lower in the case of an advection- 
dominated accretion flow (ADAF). However, there is an ad- 
ditional problem that if all of the gas captured at _Rb with 
rate M out were flowing towards the SMBH, that would con- 
tradict an upper limit on the gas density in the inner re- 
gions imposed by measurements of Faraday rotation in the 
dire ction of Sgr A*, w hich imply that M in <;10~ 7 M yr ~ 1 
fe.g. lBower et al.ll2003f ). Largely to circumvent this problem, 
radiatively inefficient accretion flow (RIAF) models have 
been proposed for Sgr A*, in which very little mass avail- 
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able at large radii actually accrete s onto the black hole (e.g. 
lYuan. Quataert. fc Nara van 2003). To summarise, the ques- 
tion of gas accretion onto the SMBH in the GC remains an 
open and hotly debated issue. 

B03 also discussed the alternative possibility of the ex- 
tended X-ray emission from Sgr A* being the sum of a large 
number of point X-ray sources. However, these authors could 
not suggest good candidates for being such objects. Below, 
we revisit this hypothesis and suggest that combined emis- 
sion of a few thousand late-type main-sequence (MS) stars 
may contribute significantly to or even dominate the X-ray 
emission from the vicinity of Sgr A*. 

The GC stellar cluster is the densest concentration of 
stars in the Milky Way. There are few 10 5 Mq of stars 
within 0.25 pc of Sgr A*, with the space density of massive 
(> IOMq) stars rising inwards down to the smallest resolv- 
able projected d i stances of ~ 0.005 pc from the SMBH (see 
lAlexanderl [20051 ; | Genzel. Eisenhauer. fe Gillessenl I20I0I for 
reviews). If the space density of presumably much more nu- 
merous late-type MS stars behaves similarly, many of these 
should be rapidly rotating as a result of repeated tidal spin- 
ups cause d by close encounters with other stars and stellar 
remnants l| Alexander fc Kumar|[200ll . hereafter AK01). We 
pursue this idea further by pointing out that fast rotation in 
stars with outer convective zones is always associated with 
strong coronal activity and greatly enhanced X-ray emission. 
Hence, if there is indeed a high-density cusp of late-type MS 
stars in the GC, which could be verified by near-infrared 
(NIR) observations with future extremely large telescopes, 
their cumulative coronal radiation can be responsible for a 
significant fraction of the X-ray emission from the vicinity 
of Sgr A*. 



2 THE MODEL 

AK01 estimated typical rotational velocities that stars in a 
GC stellar cusp can acquire over their lifespans as a result of 
repeated tidal interactions. We use their results to calculate 
the rotational evolution of stars. In addition, we take into 
account a counteracting effect - magnetic braking of stellar 
rotation. 



2.1 The stellar cusp near the Galactic SMBH 

The key question and major source of uncertainty for this 
study is how many late-type MS stars are present near Sgr 
A*. Unfortunately, NIR observations cannot yet provide a 
direct answer, since Sun-like stars are too faint and probably 
too densely clustered in the Galactic nucleus to be detected 
and individually resolved by present-day 8-10 meter tele- 
scopes. 

High spatial resolution NIR observations, which cur- 
rently can reliably resolve individual stars brighter than 
K ~ 17 in the GC region, have demonstrated that the space 
density of such bright stars increases inwards as p,(r) oc r -7 , 
with 7 w 1.2, in the innermost ~ 0.25 p c, whereas the slope 
steepens to 7 w 1.75 at larger distances (| Schodel et al 1 120071 : 
the actual slope of the nuclear cluster is inferred to be yet 
steeper once the contribution of t he Galactic bulge to ste l- 
lar counts is taken into account. iGraham fc Spitler| [2009). 



Given the limit ^<17 and taking into account the extinc- 
tion of Ak ~ 3 towards the GC, this radial profile repre- 
sents the combined space density of early-type MS stars of 
mass > 5Mq (i.e. o f B and earlier types) and low-mass g i- 
ants (see Fig. 2 in iGenzel. Eisenhauer. fc Gillessenl I201C ) . 
Furthermore, a very similar radial profile describes well the 
surface brightness of diffuse NIR light in the GC region, i.e. 
the cumulative emission of un resolved stars l|Schodel et al.l 
120071 ; lYusef-Zadeh et ai1l201lf ). The diffuse light is expected 
to be dominated by medium- (<5Mq) and low- (<LMq) 
mass MS stars. Hence, one can infer that there is a central 
cusp with 7 « 1.2 in the space density of such stars (this 
issue is further discussed in i}2.1.2l below). 

Further constraints on the radial profile and com- 
position of the GC stellar cluster are provided by spec- 
troscopic NIR observations, which can distinguish early- 
type MS stars from evolved low-mass stars, but are 
currently limited to K<15.5. Such studies have demon- 
strated that the number density of B dwarfs in the 
central ~ 0.5 pc mo notonically rises towards Sgr A* 
with a sl ope 7 ~ 2.5 llBuchholz. Schodel. fc Eckartfl2009l ; 
iDo et al l 1200 SJ; iBartko et al.l I2010T ); in particular, there 
are several tens of B stars in the central arcsecond, i.e. 
within 0.05 pc of the SMBH - the famous S-star clus- 
ter. In contrast, the surface number density of red gi- 
ants is a pproximately constant within several arcseconds of 
Sgr A* |d"o et al.ll2009l ; iBuchholz. Schodel. fc Eckartll2009l ; 
IBartko et al.l 20101 ). which implies that the radial distribu- 
tion of evolved low-mass stars is not steeper than r _1 . 

These observed distributions pro bably reflect a combi- 
nation of effects fsee lAlexanderl [20051 for a detailed discus- 
sion). Two-body interactions can establish a relaxed distri- 
bution, p*(r) oc r ~( 3 / 2 +PM) j £ s t ars i n the GC region on 
a time scale shorter than the Hubble time, with the more 
massive stars being mo re concentrated, pm ~ 0.25, than 
the lighter ones, pm ^0 (jBahcall fc Wolj|l977l ). In addition, 
the spatial distribution of early-type stars can be strongly 
affected by the recent history of star formation in the GC re- 
gion, which is not well known. Finally, the observed dearth 
of evolved stars in the innermost region may be at least 
partially caused by destruction of giants' env elopes during 
collis i ons with stars and stellar remnan ts (e.g. IGenzel et al.l 
1 19961 ; lAlexanderl H999I ; iDale et alj|2009h . 

We conclude that existing observations neither strongly 
support nor reject the possibility that there is a central cusp 
with 7 > 1 in the space density of low-mass (~ IMq) MS 
stars. 



2.1.1 Parametrisation 

Based on the above di sc ussion and following 
IGenzel. Eisenhauer. fc Gillessenl (|2010l ). we assume the 
following radial profile of total stellar mass density within 
0.25 pc of Sgr A*: 



A 



0.25 pc 



M pc" 



(1) 



where the coefficient A is discussed below. We further as- 
sume that the adopted power-law slope 7 = 1.3 pertains to 
the entire MS population in the central 0.25 pc. For com- 



pleteness, we also adopt that p* 



at r > 0.25 pc 
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IjGenzel. Eisenhauer. fc Gillessedl2010h . but this is of little 
importance for the present study, which focuses on the im- 
mediate vicinity of Sgr A*. 

The current best estimate of the coefficient of the mass 
density profile (equation [T]) is A ~ 1.5 x 1O 6 M0 pc -3 , with 
a combine d statistical and systematic uncertain ty of a fac- 
tor of ~ 2 (jGenzel. Eisenhauer. fc Gilles scn 2010). This esti- 
mate is based on dynamical measurements (proper motions 
and radial velocities of stars) of the tot al enclosed mass of 
the stellar cluster in the central parsec (jTrippe et al.ll2008l ; 
ISchodel. Merritt. fc Eckartl [2009) and the assumption that 
total mass density follows the number density of bright stars 
at r < 1 pc. It is important to emphasise that there are no 
dynamical measurements of the enclosed mass of the stellar 
cluster at distances r<0.5 pc from Sgr A*, where the gravi- 
tational potential is dominated by the SMBH, apart from a 
fairly weak upper l imit of ~ 10 5 Mq on any extended mass 
within ~ 0.005 pc (|Gillessen et al J 120091 ) . 

We further assume a ratio /lms = 0.5 for the number 
of MS stars with masses between 0.4 and 1.5 Mq (i.e. early 
M to late A stars) per Mq of total mass in the GC stel- 
lar cusp, and a ratio /sr = 0.25 for the number of stellar 
remnants (white dwarfs, neutron stars and black holes) per 
Mq of total mass. These fiducial numbers are very close to 
those adopted by AK01, who assumed a Salpeter power-law 
(a = 2.35) present-day mass function (PMF) with a low- 
mass cutoff of ~ O.4A/0. This is a good approximation for 
a continuously forming stellar population with the canoni- 
cal (such as found in the solar neighbourhood) initial mass 
function (IMF), somewhat affected by mass segregation in 
the potential well of the SMBH. As disc ussed in detail by 
lLockmann. Baumgardt. fc Kroupal (|2010l ). this scenario is 
consistent with the K-band luminosity function measured 
at K<1 7 in the innermost arcsec and also at > 12" from 
Sgr A* ijBartko et al.ll201fj ). as well as with the ratio of total 
mass to diffuse light in the central parsec. However, we stress 
that there still remains much uncertainty regarding the star 
formation history in the GC region and consequently the 
PMF of stars and stellar remnants near Sgr A*. In particu- 
lar, the disc(s) of massive Wolf-Rayet/O stars found between 
1" and 12" from Sgr A* demonstrate that a starburst with 
a highly unusual, top-heavy IMF occur ed 6 ± 2 Myr ago 
jPaumard et aLlbOOd - lBartko et alj|20ld ). 



2.1.2 Constraint from diffuse NIR light measurements 

Although, as we stressed above, Sun-like stars cannot yet be 
directly observed and counted in the GC region, it is possible 
to estimate or at least place upper limits on the number den- 
sity of low-mass stars near Sgr A* by carefully subtracting 
the contribution of individual detected stars from NIR im- 
ages and measuring the surface brightness of the remaining 
diffuse light, presumable made up by weaker, undetectable 
stars. 

Adaptive optics NIR imaging observations performed at 
the ESO VLT ijSchodel et al.ll2007l ) suggest that the cumula- 
tive surface brightness of unresolved, K ^> 17 stars increases 
approximately as IT 2 in the central few arcseconds, where 
6 is the projected distance from Sgr A* in arcsec. The mea- 
sured enclosed flux density of diffuse K-band (~ 2.3 )im) 
light is 



F„(< 9) » 0.056I 1 ' 8 Jy. (2) 

Depending on the (unknown) PMF of stars in the cen- 
tral region, this cumulative flux can be dominated by unre- 
solved stars of higher or lower mass. Conservatively assum- 
ing that at most half of the total diffuse light is provided 
by low-mass (~ 0.4-1. 5M ) MS stars, with K ~ 21-23, the 
above value of the diffuse light flux implies that there are 
<3 x 10 4 6» 18 Sun-like stars near Sgr A*, or equivalently their 
total number within 0.1 pc, 

iV L M S (r < 0.1 pc)^6 x 10 4 . (3) 

For our assumed stellar density profile (equation [TJ and the 
adopted value /lms = 0.5, this corresponds to the following 
upper limit: 

^<5 x 10 6 M© pc" 3 . (4) 

Recently. lYusef-Zadeh et alj ([201 1) studied the diffuse 
NIR light around Sgr A* using HST observations. They 
found that the surface brightness at 1.45-1.9 fim rises to- 
wards Sgr A* as 0-°-34±o.o4 and 0-o.i3io.O4 &tQ> Q 7 t, and 

8 < 0.7", respectively, and concluded that the diffuse light 
is likely dominated by low- mass (~ 0.8Mq) MS stars with a 
total number of ~ 4 x 10 4 within r = 0.1 pc. This estimate is 
in good agreement with the upper limit (equation[3]l derived 
above from the VLT AO imaging data. 

2.2 X-ray emission from spun-up stars 

2.2.1 Tidal spin-up 

AK01 presented a formalism for estimating the tidal spin-up 
of a target star of mass M» and radius R* due to a fly-by of 
another star or stellar remnant of mass m with relative ve- 
locity t>rci and impact parameter R. Both the target star and 
the impactor were assumed to be in Keplerian orbits around 
a SMBH. In addition, AK01 performed numerical SPH sim- 
ulations of close encounters, R< few i?*, when the linear 
formalism becomes inaccurate. In view of the current un- 
certainties associated with the structure and composition of 
the GC stellar cusp and with the efficiency of tidal spin-up, 
we take a simplistic approach that consists of parametrising 
our model based on the results of AK01. 

Our assumptions are as follows. First, for realistic pa- 
rameters of the GC stellar cusp, neither MS stars more mas- 
sive than few Mq nor low-mass giants, with their large mo- 
ments of inertia, can be significantly spun up by repeated 
tidal interactions. Furthemore, coronal activity driven by 
magnetic dynamo is limited to sta rs with M,<1.5Mq (i.e. 
up to late A stars, iGiidell [2004J). It is however unclear 
whether it is sub-solar mass or 1-1. 5Mq stars that will dom- 
inate in the combined X-ray emission of the stellar cusp. On 
one hand, the X-ray luminosity is a nearly constant fraction 
of the bolometric luminosity of late-type stars that rotate 
in the saturation regime (see equation 1121 below), so that 
Lx oc Lboi oc M 4 for such stars. On the other hand, the 
efficiency of tidal spin-up decreases while the efficiency of 
magnetic braking increases with increasing M* (see H2.2.2I 
below). In view of the uncertainties, we made computations 
assuming that all stars that can be spun up are Sun-like 
(M, = Mq, i?» = Rq and L* = Lq) and their space den- 
sity is /lmsP» = 0.5/9,. 

Since stars can be tidally spun up not only in mutual 
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encounters but also due to interactions with stellar rem- 
nants, the total number of impactors is defined by the sum 
/lms + /sr = 0.75 in our model. As concerns the param- 
eter /sr., white dwarfs, neutron stars and black holes are 
all expected to be present in the central cluster. We simply 
assume that all stellar remnants have Msr = M©, the same 
mass as the MS stars. 

We next need to define the stellar velocity field. Fol- 
lowing AK01, we assume that the GC stellar cusp has un- 
dergone two-body relaxati on. The one-dimension al velocity 
dispersion of stars is then (|Bahcall fc Woljl 19771 ) 



<T*(r) = 



GMbh 



Pm + 5/2 r 



(•») 



which is a weak function of stellar mass, with pm ranging 
from for the least massive stars to 1/4 for the most massive 
ones. Since we are here interested in late-type MS stars, we 
adopt pm = 0, so that 



<x*(r) w 370 



0.05 pc 



-1/2 



km s 



where we have adopted the SMBH mas s M B h = 4 x 1O 6 M 
|Ghez et al.ll200Sl ; iGillessen et alj|2008h . The distribution of 
relative velocities w re i of interacting pairs of stars (or stars 
and stellar remnants) is approximately the Maxwellian one 
with a dispersion of \/6o"* (e.g. AK01). 

We now address the question of what kind of encounters 
are important. As is detailed in AK01, the efficiency of tidal 
spin-up falls rapidly when the periseparation R p between the 
impactor and the target star becomes larger than ~ 2-3 R* . 
The spin-up resulting from a closer fly-by depends on the 
impactor's orbit in the gravitational potential of the target 
star and on the mass ratio of the two components. Crudely, 
i) the efficiency of tidal spin-up increases with decreasing 
periseparation down to R p ~ i?* when loss of mass during 
the interaction becomes important, and ii) for given R p ~ 
2R t , the net change Afi in the rotational angular velocity 
of the target star is inversely proportional to v p , the relative 
speed of the components at periastron. 

In line with our simplistic approach, we adopt that 



AD. 



An s (n h R Q /v p ) 
o 



if 
if 



R P ^ Rs 
R p > Rs 



with AD B = 0.25f2 b and Rs = 2.5R®. Here fi b = 
(GMt/Zi,) 1 ' 2 is the star's centrifugal breakup speed, which 
for solar-type stars corresponds to an equatorial speed of 
« 440 km s _1 , or an orbital period of « 3 hours (w 200 times 
the rotation rate of the Sun). We further assume that re- 
peated spin-ups Ail. have random orientations with respect 
to the target star's rotation axis ft (this, by the way, implies 
that stars occasionally can experience tidal spin-downs). 

The periastron speed v p , which enters equation (0), can 
be derived from the relative velocity of the impactor and the 
target star at infinity via the standard relation for hyperbolic 
orbits, 



Rs 



(8) 



(here it is assumed that R p = Rs). Therefore, v p ^ 
2D^(Rq I R s )Rq = 0.8f2b-R©, with the minimum reached at 
Vrel = 0, and given equation ([7]), AQ never exceeds 0.3125f2b 



in our model. In fact, since typical relative velocities in the 
stellar cusp (see equation [6| v TC \ = y/6a* > D^Rq, we only 
rarely approach this limit and typically have Afl<c0.1Slb- 

The rate of encounters experienced by a star located at 
distance r from the SMBH and lead ing to its spin-up can be 
found as (|Binnev fc Tremamdli"987l . AK01) 



dN s 
dt 



r) = 40F(/ LM s + fsK)p»(r)cr»(r)Rl 



1 + 



<rf(r) 



Re 
Rs 



(9) 



(6) djv d 



where p*(r) and <r» (r) are given by equations JT} and ©, 
respectively. Equation Q includes the gravitational focusing 
term. 

Sufficiently close to the SMBH, also collisions of stars 
leading to their destruction may become important. The 
characteristic periseparation for such catastrophic events is 
Rd ~ 0.5i?» = 0.5-R© (AK01), and the corresponding rate is 

dAd, 
dt 

QbR® \ Rq 

R^ 



4^(/lms + /sR)p.(r)a«(r)i?d 



1 + 



<ri{r) 



(10) 



Therefore, for our adopted parameter values destructive col- 
lisions occur (i? s /i?d) 2 <25 (the exact value weakly depends 
on the distance from the SMBH through the gravitional fo- 
cusing term) times less frequently than spin-up encounters. 

Our choice of values for the parameters characterising 
tidal spin-up is further discussed in H3.1.1I where we com- 
pare the results of our test simulation with the results of a 
simulation with a similar set-up by AK01. 



2.2.2 Magnetic braking 

A single spinning star cannot sustain its rotation indefi- 
nitely, because it loses angular momentum with the stellar 
wind via magnetic braking. F ollowing usual practice (see 
§13 in lSchriiver fc Zwaanl 120001 ). we parametrise this effect 
as follows: 



dQ 
lit 



mb 



if Q > D c 



(ii) 



Here, Q. c is the angular velocity at which coronal activity 
reaches saturation (see i]2.2.3l belowV We adopt fi c = O.lfib, 
which is approximately equivalent to the condition Ro' w 0.1 
for the modified Rossby number (which is approximately the 
ratio of the stellar period and the characteristic convective 
turnover time) often u sed as an act ivity saturation criterion 
for late-type stars (see lGiid"ell2004l for review). In the stan- 
dard situation of a star-forming region, the characteristic 
time tmb in equation (lll|l has the meaning of duration of 
stellar "youth" , a period when stellar activity is at its max- 
imum. Observations of stellar clusters of different ages and 
stars in the field suggest that t m b ~ 50-100 Myr for solar- 
type stars and that this time increases with d ecreasing M* , 
perhaps up to ~ 200 Myr for K/M dwarfs l|Giidell [20o3 V 
At age t 3> t mh, the decelarat ion of a star follows the Sku- 
manich law (|Skumanichlll972l ), f2(t) oc t 1 ' 2 , in agreement 
with equation (|11|) . 

Since the efficiency of magnetic braking is determined 
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by stellar mass and rotation rate, the empirical relations 
for dSl/dt pertaining to stars in the solar neighbourhood 
should be equally applicable to the GC stellar cusp, with 
the only difference being that slow-down intervals, governed 
by equation can now alternate with abrupt tidal spin- 
ups. 



2.2.3 Rotation-luminosity relation 

The X-ray luminosity of a coronall y active star depends on 
its rotation rate approximately as (Giidcl 2004) 



Table 1. Model parameters. 



L x 
Lboi 



10 
10 



if n fi c 

if fi < fi c 



(12) 



Most of the emission is at energies below 1 keV. How- 
ever, in the GC region, with its large extinction, we are more 
interested in hard X-ray emission at energies above 2 keV. 
There is a well-known trend of spectral har dening with in- 
creasi ng coronal X-ray luminosity. Following ISazonov et al.l 
( 2006) and in view of equation ([12]) . we adopt for the lumi- 
nosity of solar-type stars in the 2-10 keV energy band 



Li 



0.25 



( L * V 



(13) 



It is important to note that the above formulae, used in 
our computations, are, strictly speaking, only appropriate 
for the "quiescent" state of coronally active stars. In real- 
ity, as we discuss in detail below (in H4.2.1|l , coronally active 
stars occasionally experience giant flares during which their 
X-ray luminosity increases by several orders of magnitude 
and can become comparable to the bolometric luminosity 
of the star; moreover, the emission becomes harder. If sev- 
eral such giant flares occur simultaneously in the GC stellar 
cluster at any given time, the cumulative X-ray (2-10 keV) 
luminosity of the entire cluster can be larger by tens of per 
cent or even by a factor of few than would be expected based 
on equations (|12[) and (|13[) . 



2.3 Contribution from active binaries 

Almost certainly, rapidly rotating single stars are not the 
only kind of X-ray sources present in the GC stellar cusp. 
Among other plausible contributors to the cumulative X- 
ray emission are coronally active stellar binaries (ABs) of 
RS CVn, BY Dra and Algol types. In ABs, we deal with a 
star (in fact two stars) that rapidly rotates due to synchro- 
nisation with the orbital motion in a close binary. 

We have recently demonstrated that ABs together 
with cataclysmic variables (CVs) produce the bulk of 
the Galactic ridge X-ray emission (GRXE), the appar- 
ently diffuse background observ ed all over the Milky 
Way including the GC region ( Revnivtsev et al.| 20061; 



Revnivtsev. Vikhlinin. fc Sazonovl 



20071 ; I Revnivtsev et al.l 



2009). The combined emissivity of ABs and CVs per unit 
stellar mass in the 2-10 keV energy band is ~ 3 x 
10 27 erg s" 1 Mg\ of which ~ 1.5 x 10 27 e rg s" 1 M" 1 is 
provi ded by ABs with Lhx < 10 31 erg s" 1 (|Sazonov et al.l 
2006). Since the space density of CVs and ABs with Lh x > 
10 31 erg s" 1 is only ~ lO^M" 1 l|Sazonov et alj|2006ft . we 
expect at most a few such systems to be present within 
0.25 pc of the SMBH and none within the innermost 0.05 pc, 



Symbol, 


Explanation 


Adopted value or 


acronym 




allowed range 


A 
l \ 


IVlass density coefficient 


<.0 X 1U 1V1Q pc 


t 

JLMS 


Fraction of low-mass stars 


U.o per iViQ 


fsR 


Fraction of remnants 


0.25 per M Q 


M, 


Star's mass 


1M 




Remnant's mass 


1M 


R s 


Spin-up collision radius 


2.5R© 


Rd 


Destructive collision radius 


0.5i? o 


Af2 s 


Spin-up velocity 


0.25Q b 




Saturation velocity 


o.m b 


n 


Initial velocity 


o.m b 




Magnetic braking time 


~ 50-200 Myr 


tc 


Cluster's lifetime 


10 Gyr 


e AB 


Emissivity of ABs 


<2.5 X 10 27 erg s" 1 M" 1 


M BH 


SMBH mass 


4 x 1O 6 M 


Dgc 


Sgr A* distance 


8 kpc 



given the mass profile (equation [TJ . If we assume that the 
relative numbers of ABs and solar-type stars are the same 
as in the Solar vicinity and take into account that our 
(implicitly) adopted stellar mass function has a low-mass 
cutoff of ~ 0.4Mq (see i )2.1[) . we can expect ABs with 
Lhx < 10 31 ergs" 1 to produce e A B ~ 2.5xl0 27 ergs" 1 Mg 1 
in the GC stellar cusp. The contribution of ABs can be esti- 
mated by multiplying this fractional emissivity by the total 
mass density p*. 

In reality, there might be fewer ABs per unit stel- 
lar mass near Sgr A* than in the Solar vicinity because 
binary stars can be destroyed by encounters with single 
stars in dense stellar environments. Although luminous 
(L HX ~ 10 30 -10 31 erg s _1 ) RS CVn systems have orbital 
periods ~ 1-10 days a nd separations ~ a few 10~ 2 a.e. 
l|Strassmeier et all ll993T ) , even such tight systems can be 
evaporate d on a time s cale of 10 8 -10 10 years within 0.25 pc 
of Sgr A* <|Peretsll2009h . while one should take into account 
the fact that one or both components of RS CVn binaries 
are typically in their subgiant phase of evolution, i.e. such 
systems are usually old. Therefore, our adopted fractional 
X-ray emissivity of ABs should in fact be considered an up- 
per limit on the actual conribution of such systems. 



2.4 Final set-up 

There are a few additional parameters in our model. One 
is the lifetime of the GC stellar cusp, t c . We adopt t c — 
10 10 yr and assume that stars are formed at a constant rate 
over this time span. Another one is the rotational velocity of 
newly born stars, Qo. We adopt flo = 0.1S7b, which implies 
that stars begin their lives with marginally saturated coronal 
activity. Finall y, we assume the distance to the GC to be 
D GC = 8 kpc (|Ghez et al.ll2008l ; iGillessen et alj|2009h . 

The parameters of the model are summarised in Table Q] 
The key parameters, i.e. those whose uncertainty strongly af- 
fects the predicted X-ray luminosity of the GC stellar cusp, 
are the mass density coefficient A and, to a lesser degree, the 
characteristic magnetic braking time i m b. We allowed their 
values to vary within the indicated ranges, consistent with 
existing observations of the GC region and nearby stars, re- 
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spectively. Also, as was discussed in i)2.3l the contribution 
of active binaries, 6ab, can vary from essentially zero to the 
value quoted in Table [1] The values of the other parameters 
are fixed in the model. Most of these can be considered to be 
well determined: first of all A/bh, Dqc, and flo, and to 
a lesser degree R s , and Af2 a . The lifetime of the stellar 
cusp, t c , has very little influence on the results of the simula- 
tions. The relative fractions and typical masses of low-mass 
stars and stellar remnants (/lms, /sr., M», Msr) may actu- 
ally be significantly different from our adopted values, but 
this uncertainty is largely absorbed in the uncertainty on 
the mass density coefficient A. 

We carried out computations by combining numerical 
integration over the radial mass profile (equation [TJ with 
a simple Monte-Carlo approach for individual stars within 
a given radial bin. For each star, first its age, to, is drawn 
randomly between and t c . Then the fate of the star is 
followed from —to until t = (present day). During each 
small interval dt, the star can either spin down by mag- 
netic braking or experience a tidal impact leading to its 
spin-up. In the latter case, the star's angular velocity in- 
creases as fl(t + dt) — tt(t) + Af2, where the direction of 
Af2 is drawn randomly whereas its amplitude is calculated 
according to a relative velocity v Te \ that is drawn randomly 
from the Maxwellian distribution with a dispersion equal to 
y6a*(r) (recall that tidal spin-up is most efficient for low- 
velocity encounters). At t = 0, the X-ray luminosity of the 
star is calculated according to its final angular velocity, fi(0). 



3 RESULTS 

Below, we first discuss ( £|3.1[) the computed intrinsic prop- 
erties of cusps of spun-up, X-ray luminous stars, such as the 
number of spin-up events per star, distribution of rotational 
velocities of spun-up stars and X-ray luminosity density as 
a function of distance from the SMBH. We then compare 
(JO} the computed X-ray surface brightness profiles with 
an actual X-ray image of Sgr A*, to which end we analyse 
archival Chandra data. We then present ( H3.3P approximate 
fitting relations that allow one to predict the X-ray luminos- 
ity and radial profile of a GC stellar cusp for a given mass 
density coefficient A. Finally, in H3.4I we analyse the Chan- 
dra spectrum of Sgr A* and discuss its consistency with the 
expected X-ray spectrum of an ensemble of coronally active 
stars. 



3.1 Properties of the cusp of spun-up stars 

Figure [T] shows examples of synthetic X-ray light curves of 
stars located at different distances from the SMBH. The 
parameter values are as quoted in Table [T] in particular 
A = 4 x 10 6 M Q pc -3 and t mh = 100 Myr. We see sharp spin- 
ups caused by tidal interactions and subsequent long peri- 
ods of decaying coronal activity. Outbursts have very short 
rise timef] and characteristic durations (FWHM) ~ t m h ~ 




1 The minimal possible time scale is determined by the charac- 
teristic fly-by time ~ 2R B /v Ic i ~ 1 hour, however, the adjustment 
of the stellar magnetic field and outer layers to the new rotation 
rate will probably proceed o n the magnetic dynamo cy cle time 
scale, which is ~ 1-10 years (Saar & Brandenburg 1999). 



Figure 1. Examples of X-ray light curves for stars located at 
different distances from the SMBH, for A = 4 X 1O 6 M pc -3 , 
tmb = 100 Myr and other parameter values as given in Tabic [T] 



100 Myr. Since our model essentially ignores inefficient, dis- 
tant collisions, namely those with R p > R s = 2.5Rq (see 
equation [7)l , it typically takes just one spin-up episode to 
bring a star into the regime of coronal saturation (fi ~ Q c , 
and hence Lhx ~ 10 30 erg s _1 ). The cumulative X-ray emis- 
sion of the stellar cusp is dominated by those stars that are 
currently caught shortly after the latest spike in their light 
curve. 

Figure [2] shows the radial profiles of key average prop- 
erties of the stellar cusp. The top panel shows mass density, 
p», and X-ray (2-10 keV) luminosity density, ehx, as func- 
tions of distance from Sgr A*. The dependence enx(r) can 
be approximately described by equation (fl4)l below. The 
middle panel of Fig. shows the average number of tidal 
spin-up events experienced by stars as a function of radius. 
Finally, the bottom panel shows the radial profile of angular 
rotational velocity. Specifically, we present the average value 
(fi/fib) and a characteristic value fi/fib such that stars with 
Q ^ Cl produce half of the total X-ray emission at a given 
radius; also shown is the relative fraction of such rapid ro- 
tators. 

We see that stars experience on average (N s ) ~ 70, 
~ 20, ~ 7 and ~ 1.5 spin-up events over their lifetimes at 
0.005, 0.01, 0.02 and 0.05 pc from the SMBH, respectively. 
Outside of the innermost 0.006 pc region, the average time 
between impacts is longer than the characteristic magnetic 
braking time: (At B ) ~ 0.5t c /(N s ) > t mb - As a result, the 
majority of stars have managed to spin down after their last 
tidal interaction. However, a noticeable fraction of stars had 
their last spin-up not more than a few t m h ago and it is 
these stars that produce the bulk of the X-ray emission in 
the GC stellar cusp. For example, half of the total emission 
at 0.01 pc from the SMBH is produced by 18% of stars, 
namely those having Q J> 0.055f2b, half of the emission at 
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Figure 2. Radial profiles of average stellar properties in the GC 
cusp for A = 4 X lO 6 A/0 pc -3 , i mb = 100 Myr and other param- 
eter values as given in Table [T] Top panel: Mass density (solid 
line and left axis) and 2—10 keV luminosity density (dashed line 
and right axis). Middle panel: Number of spin-ups per star. In 
the region to the left of the vertical dotted line, stars are likely 
to be destroyed by a catastrophic collision. Bottom panel: Angu- 
lar rotational velocity (left axis): its average value (fi/fi b ) (solid 
line) and a characteristic value Q/Q\> such that stars with Q ^ f2 
produce 50% of the total (excluding the contribution of ABs) 2— 
10 keV emission at a given radius (dotted line). Also shown is the 
corresponding relative fraction of such rapid rotators (dashed line 
and right axis). 

0.02 pc by 12% of stars, with Q > 0.05Sl b , and at 0.05 pc 
by 6%, also with fi £ 0.05f2 b . 

Within 0.011 pc of the SMBH, stars are likely to be de- 
stroyed by a catastrophic collision (see the middle panel of 
Fig. [5] and H2.2.1l above). This will probably lead to a signif- 
icant depletion of the stellar cusp. In addition, by providing 
gas for accretion onto the central SMBH, such collisions can 
cause strong outbursts of Sgr A* (see H4. 3. 21 below). 

3.1.1 Comparison with AK01 

We made a special computation for the case of 10-Gyr old 
stars that begin their lives with no rotation, are subsequently 
spun up by repeated tidal interactions and do not lose angu- 
lar momentum via magnetic braking, for a stellar cusp with 
A = 3 x 10 8 Mq pc~ 3 . These conditions are nearly identical 
to the scenario considered by AKOT 

The stars acquired high rotational velocities over their 
lifespans, (fi) « 0.26fi b , « 0.17fib and « 0.10fi b at 0.02, 
0.05 and 0.1 pc from the SMBH, respectively. These val- 
ues are close to but somewhat lower than the correspond- 
ing numbers in AK01 (see their Fig. 7), ~ 0.3f2 b , ~ 0.2f2 b 
and ~ 0.15f2 b , respectively. This indicates that our choice 
of values for the R s and AQ S parameters, which define the 
efficiency of tidal spin- up in our simplified treatment, is rea- 
sonable and not overly optimistic. 




Figure 3. Chandra image (north is up and east is left) of the 
central ~ 20" X 20" of the Galaxy in the 2-8 keV energy band. 
The angular size of each pixel is 0.492" X 0.492". The positions of 
the radiosource Sgr A* and the star cluster IRS 13 are marked. 
Flares of Sgr A* have been filtered out. The circle with a radius of 
3 pixels drawn around Sgr A* denotes the region used for X-ray 
spectral analysis. 

3.2 X-ray surface brightness profile 

We now address the X-ray surface brightness profile of Sgr 
A*. Since the observations reported by B03, the Sgr A* field 
has been re-observed with Chandra many times. The total 
effective exposure accumulated is ~ 1.2 Ms, which repre- 
sents a factor of ~ 25 increase with respect to B03. We used 
these archival data to measure the X-ray surface brightness 
around Sgr A*. 

Since apart from the quiescent emission studied in this 
pape r, Sgr A* also exhibits X-ray flares l|Porquet et alj 
2008), which likely originate in the inner regions of an 
accretion flow onto the SMBH (as suggested by multi- 
wavel ength and in particular radio interferometric observa- 
tions, lYusef-Zadeh et al.| [2009 ; see however our discussion in 
i|4.2.1|l . we made an attempt to filter out such flares from 
the Chandra data. Specifically, we excluded observations in 
which increases by a factor of 4 or more in the X-ray flux 
from Sgr A* integrated over 100 s intervals were detected. 
The filtered data set has a total effective exposure of 620 ks. 
The resulting image of the central ~ 20" X 20" region of 
the Galaxy obtained with the ACIS-I detector is shown in 
Fig. El 

Even though we display only the innermost part of 
the GC region in Fig. [3j the X-ray image still exhibits a 
lot of detail apart from the source at Sgr A*. In partic- 
ular, there are two bright features within 5" of Sgr A*, 
including the complex IRS 13 of ho t, massive stars (see 
iGenzel. Eisenhauer. fc Gillessenl |2010| for a review). Both 
of these features are located to the south of Sgr A*. We 
therefore constructed two surface brightness profiles by in- 
tegrating Chandra counts in one pixel-wide semi-annuli to 
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Figure 4. Radial surface brightness profiles of the X-ray emis- 
sion to the north and to the south of Sgr A* (open and filled 
circles, respectively), in comparison with the radial profile of the 
nearby source CXOGC J174540.9-290014 (crosses) normalised 
to match the central peak of the northern profile of Sgr A*. Also 
shown are the fits of the northern profile of Sgr A* and the pro- 
file of CXOGC J174540.9-290014 by a sum of a constant and a 
Gaussian with a = 1.45 and 0.75 pixels, respectively (solid lines). 
The Gaussian components are shown separately with the dashed 
lines. 



the north and to the south of Sgr A* (Fig. |4| . For compari- 
son we show the radial profile of the presumably point-like, 
bright source CXOGC J174540. 9-290014 located ~ 20" 
from Sgr A* (and hence not shown in Fig. [3J, which was 
also discussed by B03. This profile can be fitted (taking into 
account the finite width of the pixels) reasonably well by 
a sum of a constant and a Gaussian with a w 0.37". The 
latter can be used as an approximation of the point-spread 
function (PSF) for our mosaic Chandra image of the Sgr A* 
region. Fitting the northern radial profile of Sgr A* with 
the same model gives a ~ 0.7", whereas the southern pro- 
file is significantly non-Gaussian. We can therefore estimate 
the intrinsic (i.e. PSF subtracted) size of the central source 
at w 0.6". Our estimates of the PSF width and the intrinsic 
size of the emission centred at Sgr A* are in good agreement 
with the original results of B03. 

Figure [5] compares the Chandra radial profile measured 
to the north of Sgr A* with the X-ray surface brightness pro- 
file expected for a stellar cusp with A = 4 x 1O 6 M0 pc -3 , 
tmb = 100 Myr and other parameter values as given in Ta- 
ble [1] i.e. the same model as in Fig. [2] We have converted 
all modelled profiles to Chandra counts using the spectral 
model described below in i)3.4l The dotted line is the model 
without taking into account the angular resolution of the 
telescope, while the solid line shows the model convolved 
with a Gaussian with a — 0.37" to approximate the Chandra 
PSF. We see that the unconvolved model surface brightness 
in the innermost region is approximately inversely propor- 
tional to the projected distance (see equation [15] below) . 
Finally, the dashed line shows the same model but assum- 
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Figure 5. Surface brightness profile of the extended X-ray emis- 
sion around Sgr A*. The data points are Chandra 2-8 keV mea- 
surements to the north of Sgr A*. The dotted line is a model with 
A = 4x 10 6 Mq pc~ 3 , t m t, = 100 Myr and other parameter values 
as given in Table [l] There is a ~ 20% systematic uncertainty in 
the normalisation of the model, associated with the conversion of 
the intrinsic (unabsorbed) flux to Chandra counts. The solid line 
is the same model convolved with a Gaussian with a = 0.37" , ap- 
proximating the Chandra PSF. The short-dashed line is the model 
assuming that there are no stars in the central "destruction" zone 
(see Fig. [2). The magenta dot- long-dashed line illustrates what 
would be if the stars experienced no magnetic braking. The blue 
dot-short-dashed line is a modification of the model, in which all 
stars are born simultaneously 10 Gyr ago, so that there is no 
contribution from young (and therefore coronally active) stars. 
Finally, the maximal expected contribution of active binaries is 
shown with the red long-dashed line. 

ing that the central "destruction" zone (see Fig. is devoid 
of stars. Since complete depletion of the innermost cusp by 
collisions is unlikely, the solid and dashed curves may be 
regarded as an upper and lower limits on the expected emis- 
sion, respectively. 

The model reproduces the innermost part of the Chan- 
dra surface brightness profile fairly well, underpredicting the 
measured X-ray flux from the central arcsecond by ~ 25%, 
which is comparable with the uncertainty in the conversion 
of the intrinsic (i.e. unabsorbed) source flux to Chandra 
counts (see Table [2}. For comparison, we show with the ma- 
genta dot-long-dashed line in Fig. [5] the same model, but 
assuming that the stars experience no magnetic braking. 
The resulting X-ray emission is more than an order of mag- 
nitude stronger than is observed, because all stars now spin 
with fi > O.lfit,, as a result of the rotation received at birth 
and in tidal interactions, and produce coronal emisssion with 
maximal possible efficiency. 

The effect of X-ray emission induced by tidal spin-ups 
is limited to ~ 1.5" (~ 0.06 pc) from Sgr A*. However, 
our model is also able to account for a significant fraction 
of the apparently diffuse X-ray emission detected by Chan- 
dra further away from Sgr A*. In our model, which assumes 
continuous star formation over t c = 10 Gyr, this additional 
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Figure 6. As Fig. [5] but for a set of different models. The dotted, 
solid and dashed lines are for t m b = 100 Myr and A = 2, 4 and 
6 X 1O 6 M0 pc~ 3 , respectively (collisional destruction of stars in 
the vicinity of the SMBH is not taken into account). The dash- 
dotted line is for t mb = 200 Myr and A = 3 X 1O 6 M pc -3 . 



emission is partly associated with late-type MS stars that 
were born at most several magnetic braking times (i.e. a few 
10 s yr) ago and hence continue to rotate fast and produce 
strong X-ray emission. Note that the region between 1" and 
12" from Sgr A*, where a starburst apparently occured just 
~ 6 Myr ago (|Paumard et al.ll2006l ). is expected to contain 
a large number of young, strongly X-ray emitting stellar ob- 
jects that may account for a s ignificant fraction of the appa r- 
ently diffuse X-ray emission (Navaksh in fc Sunvaevll2005l '). 

An additional significant contribution can be provided 
by active binaries (see the red long-dashed line line in Fig. [5]) 
if such systems are present and have not been evaporated 
in the GC central cusp (see the discussion in H2.3p . This 
combination of young stars and ABs resembles the situation 
in the solar neighbourhood, where young single stars, ABs 
and catacalysmic variables all pro vide similar contribu tions 
to the cumulative X-ray emission |Sazonov et al.ll2006f ). 

Importantly, the X-ray flux from the innermost ~ 1.5" 
region is essentially determined by the total mass of late- 
type MS stars in the central cusp and is almost insensitive to 
the history of star formation preceeding the present epoch. 
This is illustrated in Fig. [5] with the blue dot-short-dashed 
line, which is the result of a computation for the same model 
as before but assuming that all stars were born at the same 
time 10 Gyr ago. The flux from the vicinity of Sgr A* is 
nearly the same as in the case of continuous star formation, 
but the emission outside ~ 1.5" from Sgr A* has decreased 
due to the disappearance of young stars. 

The uncertainties associated with some of our model 
parameters (Table [TJ, in particular with the mass density 
coefficient A and magnetic braking time scale i m b, translate 
into an uncertainty of the predicted X-ray emission from 
stars in the vicinity of Sgr A*. In Fig. [5] we demonstrate how 
the X-ray surface brightness profile is modified by changing 
these parameters within reasonable ranges. We see that with 



t m b fixed at 100 Myr, the data favour (if we were to explain 
the entire X-ray emission from Sgr A* by coronal radiation 
from spun-up stars) a rather high stellar density A ~ (4- 
6) x 10 6 Mq pc -3 (taking into account the ~ 20% uncertainty 
in the model normalisation due to the conversion to Chandra 
counts). On the other hand, taking t mb = 200 Myr, which is 
still a reasonable value for <Mq stars, we can reproduce the 
emission from Sgr A* already with A ~ 3 x 1O 6 M0 pc~ 3 . 



3.3 Approximate scaling formulae 

Based on the results of our simulations, we can write down 
an approximate formula (which provides good fits to the 
computational results at least up to A — 1.2 x 1O 7 M0 pc -3 ) 
for the expected radial profile of the intrinsic, i.e. unab- 
sorbed, X-ray (2-10 keV) luminosity density of a stellar cusp 
of density p*(r) = A(r/0.25 pc) -1 ' 3 , assuming a magnetic 
braking time t mb = 100 Myr: 



eux(r) w 7.5 x 10 34 ^' 6 



+3 x 10 34 A 6 



0.05 pc 



0.05 pc 



-1 -3 

erg s pc , 



(14) 



where A 6 = A/(10 6 M© pc" 3 ). 

The first term on the right-hand side of equation (|14[) 
reflects the contribution of spun-up stars, whereas the sec- 
ond term is due to the combined emission of recently born 
(less than a few t mb ago) and hence rapidly rotating stars 
and ABs. The latter contibution obviously scales linearly 
with the cusp density. In contrast, the cumulative X-ray lu- 
minosity of spun-up stars grows approximately as A 6 with 
increasing cusp density, since it is not only the total num- 
ber of stars that increases, but also the relative fraction of 
fast rotators among them. The amplitude of the second term 
should be regarded as indicative only, because it can vary by 
a large factor depending on the actual star formation history 
in the cusp and the actual presence of ABs there. 

Given the GC distance of 8 kpc and taking into ac- 
count the absorption column Nh ~ 1.1 x 10 23 cm -2 towards 
Sgr A* (see Table [2] below), equation (|14[) translates into an 
approximate relation for the expected surface brightness of 
a GC stellar cusp in the 2-10 keV energy band: 



.(*) 



1.13 x W- ls Al- 6 6 



+8.6 x W~ 10 A 6 e~ 



erg s 



-1 -2 -2 

cm arcsec , 



(15) 



where 6 is the projected distance from Sgr A* in arcsec. For 
Ae > 2, the first term dominates at 8<l-2" from Sgr A*. 

It is important to emphasise that if magnetic braking, 
for some reason, operates less efficiently in spun-up stars in 
the GC region than in stars in the solar neighbourhood, i.e. 
tmb > 100 Myr, then rapidly rotating stars will constitute 
a larger fraction of all stars in the cusp and thus produce 
more X-ray emission per unit stellar mass, e.g. by ~ 50% if 
t m b = 200 Myr for a given value of A. A similar effect may 
result if the X-ray luminosity of the stellar cusp is domi- 
nated by (presumably) much more numerous K/M dwarfs 
rather than by Sun-like stars. In either case, it is possible 
that a somewhat lighter stellar cluster, than predicted by 
equation (|15[) , could produce a given X-ray flux from Sgr A*. 
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3.4 X-ray spectrum 

We have demonstrated that a high-density cusp of late-type 
MS stars in the GC region can produce X-ray emission with 
the radial surface brightness profile as observed by Chandra 
near Sgr A*. Will this emission from rapidly spinning stars 
also have a spectrum similar to that of Sgr A*? 

To answer this question, we re-measured the spectrum 
of Sgr A* reported by B03. Since the GC region is rich 
in high-energy phenomena, a significant or even dominant 
fraction of the X-ray emission observed outside of the in- 
nermost few arcseconds (see the image in Fig. [3]) may have 
a different nature from the central X-ray cusp. We there- 
fore constructed a spectrum from the photons detected 
within 3 pixels (w 1.5") of Sgr A*. We modelled the de- 
tector background for Chandra/ ACIS using the acisbg task 



I http://cxc.cfa.harvard.edu/contrib/maxim/acisbg/ ) based 
on the stowed dataset. The background-subtracted 2-8 keV 
flux from the studied region is 1.45 x 10 -13 erg s _1 cm -1 
(with a statistical uncertainty of a few per cent), in good 
agreement with B03. 

In agreement with B03, the measured spectrum of Sgr 
A* (Fig. [7J is consistent with thermal emission of optically 
thin, hot plasma, modified at low energies by absorption in 
cold gas. In addition, we detect a ~ 3<j significant excess 
that is consisent with a 6.4 keV fluorescent line of neutral 
iron. We fitted the Chandra data in the 2-8 ke V energy 
band by different models in XSPEC (version 12.5.1. lArnaudl 
1996) and found that a good fit is provided by the model 
wabs * (cevmkl + gauss). The best-fit parameters are given 
in Tableland the sp ectral fit is shown in Fig. [7] CEVMKL 
is a s um of MEKAL llMewe. Groncnsc hild. fc van den Oordl 
ll985l ; lLiedahl. Osterheld. fc Goldsteinlll995l ) models for op- 
tically thin plasmas with a range of temperatures and 
a power-law distribution of emission measures, dEM oc 
(T/r max ) Q " 1 dT/T max , where we fix a = 1 (the results 
are only weakly sensitive to this parameter). In fact, using 
a single-temperature MEKAL model instead of CEVMKL 
provides a similarly good fit to the strongly absorbed spec- 
trum of Sgr A*, but we have chosen CEVMKL because it is 
more physically justified in application to X-ray spectra of 
coronally active stars. 

Table [2] also gives the absorption-corrected flux and lu- 
minosity of Sgr A* in the 2-10 keV band for the adopted 
spectral model. These values have ~ 20% uncertainties. 
In summary, all of the characteristics of the X-ray spec- 
trum of Sgr A* derived here are in good agreement with 
the estimates of B03 (although these authors used a single- 
temperature thermal model), but the associated statistical 
uncertainties have now become small due to the greatly in- 
creased observational time and statistics. 



3.4-1 Comparison with stellar spectra 



According to the reference mass density model (equation [TJ 
A = 4 x 1O 6 M pc -3 ), which we used to describe the X-ray 
surface brightness profile (Fig. O, the total mass of stars 
and stellar remnants within the cylinder with a projected 
radius of 1.5" around Sgr A* is 8.4 x 1O 4 M (6.3 x 1O 4 M ) 
within ±0.25 pc (±0.1 pc) along the line of sight. Since 
within ~ 0.1 pc from Sgr A* the bulk of the coronal X- 
ray luminosity is produced by ~ 10% of low-mass stars (see 
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Figure 7. Spectrum of X-ray emission from the region with 
radius of 1.5" around Sgr A* measured by Chandra in the 
2-8 keV energy band (data points) and its best fit by the 
wabs * (cevmkl + gauss) model with the parameters given in Ta- 
ble [2] For comparison shown are the best- fit cevmkl models (see 
Table [2} of ASCA spectra of the coronally active stars V711 Tau 
(blue dotted) and 47 Cas B (red dashed), modified by cold gas ab- 
sorption with iVjj = 1-08 x 10 23 cm -2 and normalised to match 
the spectral flux of Sgr A* at 3 keV. The inset shows a blow- 
up of the Sgr A* spectrum at ~6-7 keV. The 6.4 keV line and 
CEVMKL components are shown by the dotted and dashed lines, 
respectively. 



the bottom panel of Fig. [5] and the related discussion in 
and we assumed /lms = 0.5, the X-ray flux accumulated by 
Chandra within 1.5" of Sgr A* can be produced by some 
3-5 thousand rapidly spinning stars, each emitting 10 29 - 
10 30 erg s _1 in the 2-10 keV energy band. 

As was mentioned in £12.2.31 there is a well-known trend 
of X-ray spectral hardening with increasing coronal lumi- 
nosity, which in turn is closely linked to stellar rotation. 
Although the rotational velocities are expected to vary from 
one star to another in a GC stellar cusp, its cumulative X- 
ray emission will be dominated by stars that are in (fi > fi c ) 
or close to being in (f2<jO c ) saturation regime of coronal ac- 
tivity (see again the bottom panel of Fig. 2) . 

An ultimate spectral test of our model would consist in 
constucting a composite X-ray spectrum of an ensemble of 
stars rotating at various angular velocities representing the 
distribution of Q values predicted by the model. We feel that 
performing such a computation is premature given the rela- 
tive simplicity of our current model. In particular, this model 
deals with only one type of star, namely IMq dwarfs. In real- 
ity, depending on the composition of the GC stellar cusp, its 
X-ray luminosity may be dominated by K/M dwarfs rather 
than by less numerous Sun-like stars. The rotational velocity 
distributions and X-ray spectra may be somewhat different 
for stars of different types, and a composite spectrum of the 
stellar cusp should take this into account. Perhaps more im- 
portantly, the combined spectrum of the stellar cusp may 
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Table 2. Best-fit spectral parameters for the X-ray emission within 1.5" of Sgr A* and for two coronally active stars. 
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Parameter 


Sgr A* 


V711 Tau 


47 Cas B 




Chandra/ACIS-I (2-8 keV) 


ASCA/GIS (2-9 keV) 




wabs * (cevmkl + gauss) 


cevmkl 


T max (keV) 


3.2 ±0.3 


3.74 ±0.08 


2.4 ±0.3 


a 


1 (fixed) 


1 (fixed) 


1 (fixed) 


A Fc 


0.80 ±0.12 


0.32 ± 0.07 


1 (fixed) 


Nn, cm" 2 


(10.8 ±0.6) x 10 22 






EW 6 . 4 kcV (eV) 


110+ 30 






X 2 /dof 


156/130 


419/315 


33/49 


Unabsorbed 2-10 keV flux (erg s — 1 cm -2 ) 


(4.3 ±0.8) x lO" 13 


2.9 x lO" 11 


2.4 x lO" 12 


Distance (pc) 


8 x 10 3 


29.0 


33.6 


Unabsorbed 2-10 keV luminosity (erg s" 1 ) 


(3.3 ±0.6) x 10 33 


2.9 x 10 30 


3.2 x 10 29 



V711 Tau is indicative only, given the simplicity of the adopted mode l. The quoted dist ances to the stars are based on Hipparcos 

parallax measurements jPerrvman fc ESAlll997T l. 



contain a significant contribution from giant stellar flares 
(with Lx ~ -t/boi), especially at energies above ~ 3 keV. As 
discussed in £|4.2.1l below. the statistics of such flares remains 
very limited. Moreover, such events are usually detected by 
all-sky monitors and their X-ray spectra are poorly studied 
at brightness peak. This renders it practically impossible to 
estimate the contribution of giant flares to the cumulative 
X-ray spectrum of the stellar cusp. In view of these compli- 
cations and uncertainties, we restrict our discussion below 
to a fairly qualitative comparison of the spectrum of the ex- 
tended emission from Sgr A* with the spectra of a couple 
of well-studied coronally active stars in the Solar neighbour- 
hood. 

The first object, V711 Tau, is a well-known RS CVn 
system. It consists of KO-2 IV an d G0-5 V s tars, tidally 
locked with a period of 2.84 days (|Fekellll983l ). Given the 
masses and radii of the components of Mi = 1.4 ± Q.2Mq, 
M 2 = 1.1 ± 0.2Mp , Ri = 3.9 ± O.2i? and R 2 = 1.3 ± 
O.27i jFekellll983l ). the subgiant component rotates at a 
significant fraction of its breakup velocity, fi/fib ~ 0.27, and 
is therefore definitely in saturated regime of coronal activity. 
We emphasise that although V711 Tau is an active binary, 
its spectral properties should resemble those of a single star 
that has just been spun up by a strong tidal interaction to 
fi > fi c . 

The other object, 47 Cas B, is a relatively young (~ 
100 Myr old) solar analog (spe ctral type GO- 2 V) be longing 
to the Pleiades moving group dTelleschi et al.ll2005l). It has 
a likely rotation period of ~ 1 day (see iTelleschi et al .1120051 
and references therein), so that fi/fib ~ 0.1, with this fast 
rotation being due to the star's young age rather than due to 
its binarity (the star is actually a member of a wide binary 
system). It might thus be close to but not necessarily fully in 
saturation regime of coronal activity. Therefore, 47 Cas B 
probably more closely resembles a typical spun-up star in 
the GC stellar cusp than V711 Tau. 

We obtained the X-ray spectra of these two systems 
using archival data of the GIS instrument aboard the ASCA 
observatory. Both spectra are well fit in the 2-9 keV band by 
the CEVMKL model. The best-fit parameters are given in 
Tabled The fit for V711 Tau could be improved further by 
allowing the abundances of S, Si and other metals to vary, 
but such small modifications are unimportant for our current 



purpose of comparing the spectral continua of Sgr A* and 
coronally active stars. Fig. [7] shows the best-fit models for 
V711 Tau and 47 Cas B, modified by cold gas absorption 
with Nn = 1.08 x 10 23 cm -2 , as measured for Sgr A*, and 
normalised to the flux of Sgr A* at 3 keV. 

We see that the spectra of V711 Tau and 47 Cas B 
are slightly harder and somewhat softer, respectively, than 
the spectrum of Sgr A*. This appears to be consistent with 
i) V711 Tau having saturated coronal activity - as wit- 
nessed by its very high hard X-ray luminosity, Lhx ~ 3 x 
10 30 erg s" 1 (2-10 keV, Table© and ii) 47 Cas B being not 
quite satura ted - its X-ray lumin osity of ~ 2.5 x 10 30 erg s _1 
(0.1-10 keV. lTelleschi et alj2005l . of which ~ 3x 10 29 erg s" 1 
is at 2-10 keV, Table© might be somewhat below the sat- 
uration level of ~ 10~ 3 Lboi- 

Although it is difficult to draw firm conclusions from 
this exemplary analysis, it indicates that the non-flaring 
emission of an ensemble of spun-up stars can produce a 
spectrum that is sufficiently hard to explain the low-energy 
(below ~ 4 keV) spectral continuum observed from Sgr A*. 
However, to account for the hard X-ray (at energies above 
~ 4 keV) flux from Sgr A* within the stellar cusp model, one 
possibly needs an additional significant contribution from 
giant stellar flares. This appears plausible (see £14.2. 1 1 and 
i ]4,2.2,3[) but requires further study. 



3.4.2 Iron abundance 

The Chandra spectrum of Sgr A*, in particular the strength 
of the 6.7-7.0 keV blend of FeXXV and FeXXVI lines, indi- 
cates a somewhat sub-solar (~ .7-0.9 relative to the solar 
value of iGrevesse fc Sauval|[l998l ) abundance of iron in the 
X-ray emitting plasma. This result appears to be consistent 
with the model of cumulative emission of coronally active 
stars. 

Indeed, both V711 Tau and 47 Cas B exhibit signif- 
icantly subsolar iron abundances. For V711 Tau, this is 
already suggested by the ASCA spectrum presented here: 
Af c ~ 0.3. More reliable estimates are provided by anal- 
yses of high-resolution spectra obtained with the Reflec- 
tion Gratin g Spectrometer on X MM-Newton, which indicate 
A Fe ~ 0.2 ijAudard et al.ll2003l ) and ~ 0.5 (|Telleschi et all 
|2005|) for V711 Tau and 47 Cas B, respectively (for the Fe 
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abundance of iGrevesse fc SauvaJll998h . These examples re- 
flect the well-known phenomenon of reduced abundance of 
iron and other low first ionisation p otential elements in ac- 
tive stellar coronae fsee lGiideill2004l for a review). 



3.4.3 6.4 keV emission 

Combined radiation of coronally active stars can also possi- 
bly account for at least part of the flux of 6.4 keV emission 
apparent in the Chandra spectrum of Sgr A* (see Fig. [7|). 
Indeed, a significant fraction (up to 50%) of the hard X- 
ray emission produced in the corona of a star irradiates the 
underlying chromosphere and photosphere and produces flu- 
orescent emission, mainly in the 6.4 keV line of iron. Such 
fluorescent emission provides a means of probing the geome- 
try of coronal flares and the photospheric abundance of iron. 

This problem has been of interest fo r a long time in 
appli cation to the Sun, first theoretically l)Ba il ll979l : iBaskol 
1 19791 ) and then observationally as powerful 6. 4 keV emission 
was actually observed dur ing solar flares (e.g. ICulhane et al.l 
Il98ll : IPhillips et alJll993 ). Recently, 6.4 keV emission has 
also been detected during coronal fl ares of few other s tars, 
including the RS CVn b inary II Peg (lOsten et al.ll2007n . the 
G-type giant HR 9024 (iTesta et al.ll2008l ) and the M dwarf 
EV Lac (Ostenetal .120101 ). The measured equivalent widths 
of the 6.4 keV line range from a few tens to ~ 200 eV, which 
is similar to the values reported for solar flares and agrees 
with the expected flux of fluorescent emission arising due to 
irradiation of a stellar atmosphere of solar Fe abundance by 
hard X rays fro m a coronal loop located at <,0.3 R t above 
the atmosphere l|Drake. Ercolano. fc Swart j|2008h . 

To our knowledge, there have been no reports of 6.4 keV 
emission in non-flaring spectra of coronally active stars, 
probably due to the insufficient sensitivity and/or low spec- 
tral resolution of the existing instruments and the presence 
of the much stronger 6.7 keV and 7.0 keV line complexes 
of iron in coronal spectra. However, based on the existing 
evidence for solar and stellar flares we may expect the inte- 
grated spectrum of a cusp of rapidly spinning stars to con- 
tain a 6.4 keV line with an equivalent width EW ~ 50 eV. 
Moreover, as is further discussed below in i]4.2.2.2l if gi- 
ant coronal flares contribute significantly to the cumulative 
X-ray emission of the cusp, somewhat stronger fluorescent 
emission can be expected, possibly consistent with the value 
EW = II0I4Q e V measured in the spectrum of Sgr A*. 

Another possibility to explain part of the observed 
6.4 keV emission is fluorescent emission from a cold gas 
permeating the stellar cusp. However, in order to account 
for, say, EW ~ 50 eV, this gas must have a column den- 
sity ~ 3 x 10 23 cm -2 , assuming irradiation by ~ 3.2 keV 
thermal emission from Sgr A*, solar abundance of iron 
(Fe/H = 3. 16 x 10~ 5 ) and adopting t he ionisation cross sec- 
tions from IVerner fc Yakovlevl (|l995T l and fluorescent yield 
of 0.34. This is much larger than the total absorption column 
of m 10 23 cm -2 measured by Chandra towards Sgr A*. 

We conclude that photospheric fluorescent emission ac- 
companying the coronal activity of spun-up stars can possi- 
bly explain the tentatively detected 6.4 keV emission from 
Sgr A*. 



4 DISCUSSION 

We have demonstrated that the combined coronal radia- 
tion of a high-density cluster of late-type MS stars "reju- 
venated" by tidal interactions can in principle explain the 
main properties of the quiescent X-ray emission from Sgr A*, 
namely its characteristic size (FWHM ~ 1-5"), luminos- 
ity (~ 10 33 erg s _1 , observed in the 2-10 keV band) and 
possibly (this issue needs further study) spectrum (opti- 
cally thin thermal emission with fcT<3.5 keV, absorbed by 
Nb. ~ 10 23 cm -2 ). In this scenario, extended X-ray emission 
is produced by several thousand rapidly spinning stars (with 
periods of > 1 day) with nearly saturated coronal activity. 

Below, we first discuss the implications of our model for 
the mass density of the GC stellar cusp ( t|4.1jl . then compare 
the observational properties of the model with those of truly 
diffuse X-ray emission from a hot accretion flow onto the 
SMBH f M4.2p and finally discuss the implications for the 
accretion onto the SMBH (JO|. 



4.1 Implications for the GC stellar cusp 

The existence of a cusp in the number density distribution 
of low-mass MS stars near Sgr A* has not yet been proven 
by direct observations. According to our baseline model, to 
explain the bulk of the X-ray emission from Sgr A* there 
must be ~ (2-3) xlO 4 Sun-like stars within a column of 
1" radius centered on Sgr A* (within 0.25 pc along the 
line of sight). This number corresponds to a range ~ (4- 
6)x10 6 Mq pc -3 for the coefficient A characterising the ra- 
dial profile of stellar mass density (equation [1} and to a 
relative fraction /lms = 0.5 of low-mass stars per Mq. 

Equations (I14|l and (|15p can be used to predict the con- 
tribution of coronal radiation of low-mass stars (not only 
spun-up stars but also young stars and active binaries) to the 
observed X-ray emission from the Sgr A* region for a given 
density of the GC stellar cusp, assuming that tidally spun-up 
stars experience efficient magnetic braking (t m b = 100 Myr). 

As follows from our discussion in EI2.1I the space 
density of stars required to explain the quiescent X-ray 
flux from Sgr A* somewhat exceeds the range A ~ (1— 
3) x10 6 Mq pc -3 , suggested by the available statistics of 
bright stars in the central parsec of the Galaxy. Further- 
more, it is possible that low-mass stars near Sgr A* are 
underabundant near Sgr A*, as suggested by the paucity of 
red giants in the central several arcseconds. However, the 
number of low-mass stars in the central arcsecond required 
to explain the Chandra data is consistent with the upper 
limit deduced from the measured surface brightness of NIR 
diffuse light near Sgr A* (see equations [3] and [4} . 

Furthermore, it is possible that a cusp of low-mass stars 
will produce a higher X-ray luminosity for a given density 
than implied by our baseline model: 

(i) First, if magnetic braking operates somewhat less effi- 
ciently than assumed, then a larger fraction of the stars will 
retain their fast rotation since their last spin-up episode. 
As we discussed in £13-31 increasing the characteristic time 
£ m b from 100 Myr to a still reasonable (especially for K/M 
dwarfs) 200 Myr leads to a 50% increase in the luminosity 
of the stellar cusp. 

(ii) Secondly, as discussed in ^4.2.11 below, multiple giant 
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stellar flares may significantly increase the X-ray luminosity 
of the stellar cusp. 

(iii) Thirdly, in our baseline model stellar remnants 
played a relatively unimportant role, as we assumed their 
relative fraction to be /sr = 0.25 and the average mass 
Msb., = IMq. In reality, the innermost region of the 
Gal axy may be overpopulated by stellar-mass black holes 
(e.g. Morris 1993; Frcitag, Amaro-Scoanc, & Kalogcra 2006; 
lLockmann. Baumgardt. fc KroupalbOlOT ). These heavy ob- 
jects can spin up low-mass stars much more efficiently than 
stars or white dwarfs (see AK01). 

(iv) Finally, in our simulations we implicitly assumed that 
stars were in circular orbits around the SMBH. A more accu- 
rate computation should take into account mutual interac- 
tions of stars moving in eccentric orbits around the SMBH, 
which might increase the number of spun-up stars. 

In summary, given the remaining uncertainties on the 
observational side and the fact that the presented model is 
somewhat simplistic, it seems plausible that a factor of few 
less stars in the central region might still produce the bulk 
of the X-ray flux observed from Sgr A*. 

4.2 Comparison with accretion flow models 

The crucial question is how to distinguish observationally 
the combined coronal radiation of a cluster of spun-up stars 
from the diffuse emission of a quasi-spherical accretion flow 
of hot gas onto the central SMBH, which has been the 
favoured explanation for the quiescent X-ray emission from 
Sgr A* so far. Below, we discuss differences and similarities 
between the two models. 

4-2.1 Variability and flares 

We have so far assumed that the X-ray luminosity emitted 
by a coronally active star is solely determined by its cur- 
rent rotation rate. In reality, virtually all stars with coronae 
experience flares during which their X-ray luminosity in- 
creases by orders of magnitude. Therefore, the cumulative 
X-ray emission of a GC stellar cusp must be variable at some 
level. 

Statistics of coronal flares i n stars of diff erent types im- 
plies (see §13 in the review of iGiidell 12004 and references 
therein) that X-ray/UV flare luminosities (or total energies) 
obey a power law with a slope of « 2, i.e. dN/dL oc L~ 2 , and 
it remains unclear if the strongest or weakest flares domi- 
nate in the long term-integrated energy output of stellar 
coronae and if the "quiescent" luminosity is in fact just a 
superposition of weak flares. 

Particularly interesting are giant flares observed 
from nearby RS CVn systems. Only a few dozen such 
flares have been detected in several decades of X- 
ray astronomical observations by all -sky monitors (see 
lArefiev. Priedhorskv. fc Borozdir] 12003). In particular, the 
SSI instrument aboard the Ariel V spacecraft detected 
7 flares that were attributed to RS CVn binaries 
l|Pve fc McHardvlll983h , including 2 flares from one of the 
best studied systems II Peg. Interestingly, all these flares 
had similar properties, namely peak luminosities of (1— 
6)xl0 32 erg s _1 (in the 2-18 keV energy band), charac- 
teristic durations ~ 5 hours and total emitted energies (in 



the X-ray band) ~ 6 x 10 36 erg s 1 . Based on this statis- 
tics, iPve fc McHardvl (|l983l ) concluded that a typical RS 
CVn binary produces about 1 such giant flare per year. This 
explains why dedicated observations of individual RS CVn 
systems by X-ray telescopes, which typically last less than 
a day, fail to detect such extreme flares. 

On Dec. 16, 2005, another giant flare was detected 
from II Peg by the hard X-ray instrument BA T and X-ray 
telesc ope (XRT) aboard the Swift spacecraft (|Osten et al.l 
I2007T ). These observations provided uniquely detailed, 
broad-band information on a giant flare. During the flare, 
the luminosity peaked at ~ 1.3 x 10 33 and ~ 8 x 10 32 erg s" 1 
in the 0.8-10 keV and 10-200 keV, respectively. Thus, the 
total X-ray luminosity of the corona for 1-2 hours reached 
;> 40% of the bolometric luminosity of the binary system 
(~ 5 x 10 33 erg s" 1 . iMewe et al.lll997l ). For comparison, 
the non-flaring X-ray (0.5-12 keV) luminosity of II Peg is 
~ 10 31 erg s _1 IjHuenemoerder. Canizares. fc SchulzH200ll ). 
or ~ 10~ 3 of its bolometric luminosity. This implies that 
stellar coronae that are already in a saturation regime can 
brighten up by an additional factor of J> 100 during giant 
flares. 

New data on giant flares are now coming from the 
MAXI all-sky monitor on the International Space Station. 
Already during the first year of operation it detected 14 
flares from 7 RS CVn binaries, including 2 huge flares with 
peak (1.5-10 k eV) luminosities ~ 5 x 10 33 erg s _1 from II Peg 
and GT Mus l|Tsuboi et alj|201lf ). These new results, first, 
confirm that giant flares occur roughly once a year per RS 
CVn system and, second, suggest that stellar coronae may 
produce even more extreme flares than was considered pos- 
sible before. 

The spun-up stars in a GC cusp should be similar to 
RS CVn systems as regards their coronal activity, apart from 
being less luminous by a factor of few, since the X-ray bright 
component in RS CVn binaries is often a subgiant rather 
than a normal star. We may thus expect every spun-up star 
to undergo a giant X-ray flare with a luminosity of 10 32 - 
10 33 erg s _1 and duration of a few hours roug hly once a 
year. A "superflare" of the M dwarf EV Lac (|Osten et al.l 
i2010h . already mentioned in H3.4.3I at the peak of which 
the X-ray (0.3-100 keV) luminosity reached more than 3 
bolometric luminosities of the star, might serve as a good 
example of such a giant flare on a single coronally active star 
(note, however, that this particular flare was fairly short, 
~ 0.5 hour). 

Given that some 3-5 thousand spun-up stars are re- 
quired to produce the measured X-ray flux from Sgr A*, we 
may expect several stars to experience giant flares at any 
moment. The cumulative luminosity of these flares may be 
comparable or exceed the total luminosity of the remaining 
several thousand "quiescent" stars. In fact, due to flaring ac- 
tivity the integrated X-ray luminosity of a GC stellar cusp 
can be significantly higher than expected based on equa- 
tions |l2l) and (fT3)> . 

If the cumulative emission of spun-up stars is domi- 
nated by strong flares, we can expect the X-ray flux from 
Sgr A* to vary on hourly and daily time scales by tens of 
per cent or even by a factor of few, due to the changing 
number and intensity of flares. This variability should be 
especially strong at higher energies, i.e. around 10 keV and 
above, because the coronal emission spectrum hardens dur- 
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ing flares, as for example wa s the case during the superflare 
of II Peg (|Osten et alj|2007r i. In principle, the already exist- 
ing Chandra and XMM-Newton data can be used to verify 
this prediction. Furthermore, it is possible that some of the 
weak er, ~ 1 hour long X - ray flares observed from Sgr A* 
fe.g.. lBaganoff et alj|2003 ; iPorquet et alj|2008| ) were in fact 
produced by stellar coronae. 

At the same time, it is clear that the stronger X-ray 
flares observed from Sgr A*, during which its luminosity 
reaches several xlO 34 and even more than 10 35 erg s _1 , are 
of different origin and most likely produced in the close vicin- 
ity of the SMBH (see lGenzel. Eisenhauer, fc Gillessenll2010l 
for review). In particular, by analysing the Chandra data 
for some of such flares we have verified that their origin is 
consistent with a point source (within the Chandra PSF) at 
the position of Sgr A*. 



4-2.2 Multifrequency properties 

The primary radiative component of active stellar coronae 
is optically thin X-ray emission of a multitemperature gas, 
with a dominant temperature kT ~ a few keV. A similar 
X-ray spectrum can be produced by a hot accretion flow 
if the accretion rate is sufficiently low. In this case, the 
bulk of the total X-ray luminosity is due to gas near the 
Bondi radius Rb ~ 0.05 pc, whereas the additional emission 
produced in the inner regions of the presumably advection- 
dominat ed accreti o n flow is relatively unimportant. As ar- 
gued bv iQuataerd |2002) , such a situation may be realised 
in the quiescent state of the SMBH associated with Sgr A*, 
for which Lhx ~ 10 _11 I/Edd, where Lhx and L^dd are the 
observed X-ray luminosity and the SMBH Eddington lumi- 
nosity, respectively. 

Given the same thermal character of the X-ray emission 
in the scenarios of a cluster of spun-up stars and hot accre- 
tion flow onto the SMBH, are there spectral features that 
could distinguish these models from one another? 

4.2.2.1 Thermal X-ray lines Both models naturally 
produce coronal X-ray lines, in particular a 6.7 keV line 
of FeXXV. As we showed in £13.4.21 the Chandra spectrum 
in the 6.7 keV region indicates a somewhat subsolar abun- 
dance of Fe for Sgr A*, which is in line with the known 
trend for active stellar coronae. However, the measured devi- 
ation from the solar abundance is minor (Af 6 ~ 0.7-0.9) and 
does not seem to contradict the hot accretion flow model for 
Sgr A* either. Indeed, recent studies of chemical abundances 
in the atmospheres of supergiants belonging to young stellar 
clusters near Sgr A* have demonstrated that these several 
million-old stars were formed from an interstellar med ium 
with a nearly solar Fe/H ratio (see iDavies et all 120091 and 
references therein). Therefore, one can expect the gas cur- 
rently feeding the SMBH to have approximately solar Fe 
abundance as well. 

What could the widths of the X-ray lines tell us? In 
the stellar cusp model, the line widths should reflect the 
emission-weighted, radially integrated distribution of pro- 
jected stellar velocities. Given equation ([6]) and the model 
surface brightness profile shown in Fig. [5] the lines pro- 
duced within ~ 1" of Sgr A* should be broadened by 
cr ~ 500 km s" 1 , which corresponds to FWHM ~ 25-30 eV 



for the 6.7 line of iron. Hence, the lines are unresolvable with 
CCD detectors. 

It is more difficult to predict the line widths in the case 
of a hot accretion flow onto the SMBH. In ADAF-type mod- 
els, ions tend to acquire the virial v elocities in the grav i- 
tational potential of the black hole (|Naravan fc Yill"l994r ). 
However, in the case of Sgr A* most of the luminosity mea- 
sured by Chandra from the central arcsecond is believed to 
be produced near the Bondi radius. In this case, the broad- 
ening of the lines of heavy elements is probably determined 
by bulk motions in the gas, which should be of the order of 
the stellar velocity dis persion if the gas is supplied by win ds 
of hot stars (see, e.g., lQuataertll2004 ; ICuadra et alJI^OOrj ). 

Therefore, to a first approximation, thermal X-ray lines 
should be similarly broadened in both scenarios. 

4.2.2.2 6.4 keV line A stand-alone issue is 6.4 keV 
emission of neutral (or weakly ionised) iron. As we discussed 
in i]3.4.3l we have tentatively detected such a line, with an 
equivalent width EW = H0I4Q e ^ in the Chandra spec- 
trum. 

There is no obvious reason to expect 6.4 keV emission 
in the scenario of a hot accretion flow onto the SMBH. On 
the contrary, as we discussed in H3.4.31 a 6.4 keV line is 
naturally produced in stellar coronae, due to the irradia- 
tion of the underlying photosphere by the coronal radiation 
and subsequent fluorscence. We noted, however, that "qui- 
escent" coronae are unlikely to produce 6.4 keV emission 
stronger than EW ~ 50 eV. However, during strong flares, 
the coronal spectrum becomes significantly harder (see be- 
low), hence a larger fraction of the total luminosity is emit- 
ted above the ionisation thereshold of iron (7.1 keV) and 
consequently a 6.4 keV line with a larger equivalent width 
can be produced. For example, during a "superflare" from 
the M dwarf EV Lac, 6.4 ke V emission with E W ~ 200 eV 
was detected by Swift/XRT (|Osten et al.ll2010l ). 

Therefore, if strong flares significantly contribute to the 
X-ray luminosity of the stellar cusp, we may expect a strong 
6.4 keV line with EW<100 eV to be present in the X-ray 
spectrum. This feature, if confirmed with higher statistical 
significance by future observations (in particular, by long 
observations with the Chandra High Energy Tranmission 
Grating, HETG), is thus a clear indicator in favour of the 
stellar cusp model for the extended X-ray emission from 
Sgr A*. 

4.2.2.3 Hard X-ray emission If the X-ray source at 
Sgr A* is associated with a Bondi/ ADAF accretion flow 
onto the SMBH, then in addition to the dominant soft X- 
ray emission originating in the vicinity of the Bondi ra- 
dius there should be a hard X-ray continuum due to the 
bremsstrahlung emission from the much hotter accretion 
flow at R <C -Rb, whose luminosity (at energies above 
10 keV) might be ~ 1/3 of the X-ray luminosity (below 
10 keV) of Sgr A* (|Quataerdl2002h . Non-thermal (e.g. syn- 
chrotron or syncrotron-self-Compton) emission produced in 
the innermost region of the accretion flow may provide some 
additional hard X-ray flux. Therefore, a hard X-ray contin- 
uum is a generic feature of Bondi/ ADAF models for Sgr A*. 

A significant hard X-ray continuum is also expected in 
the case of a cluster of spun-up, coronally active stars. In- 
deed, as we have already discussed, strong coronal flares can 
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contribute significantly to the tot al X-ray luminosity of the 
GC cusp. It is well -known (see e.g. iDennis fc Schwartdl 19891 ; 
lAschwandenll2002l for reviews) that in solar flares only part 
of the total energy released via magnetic reconnection goes 
into heating of thermal plasma, and only a fraction of this 
thermal energy is radiated away as X-ray (and UV) emission. 
A significant amount of energy goes into acceleration of non- 
thermal electrons and ions. The non-thermal electrons lose a 
small fraction of their energy as hard X-ray bremsstrahlung 
and radio gyrosynchrotron emission. 

During solar flares, non-thermal hard X-ray emission 
is some 5 orders of mag nitude weaker tha n thermal X-ray 
emission (see Fig. 94 in lAschwand cn 2002). However, dur- 
ing the 2005 giant flare of the RS CVn binary II Peg, its 
hard X-ray and soft X-ray luminosities were comparable 
(see 34.2.11 above). Therefore, if strong flares play a signifi- 
cant role in the combined emission of the GC stellar cusp, 
it might emit in total <10 33 erg s _1 in hard X-rays (at en- 
ergies above 10 keV), corresponding to a hard X-ray flux 
<ClO~ 13 erg s _1 cm -2 for the GC distance of ~ 8 kpc. 
Such a signal can in principle be detected by the JVuS- 
TAR telescope, to be launched in 2012, which will have 
~ 10" (FWHM) angular resolution and reach sensitivity 
~ 10~ 14 erg s _1 cm " 2 in the 10-30 keV energy band for 
an exposure of 1 Ms (|Harrison et al.ll2010l ). As we discussed 
in £|4.2.1l above. this hard X-ray signal should be significantly 
variable on time scales of hours and days as a result of flaring 
activity on stars. 

We conclude that both the stellar cusp and the accretion 
flow models predict hard X-ray bremsstrahlung emission at 
some significant level, although these predictions are very 
uncertain at the moment. 

4.2.2.4 Gamma-ray emission Besides electrons, also 
protons and ions are accelerated in coronal flares, reach- 
ing kinetic energies up to GEVs. In the case of the Sun, 
just a small fraction of such sub-cosmic rays escape into 
interplanetary space, while the majority eventually hit the 
lower corona and chromosphere and lose en ergy there (see 
IDennis fc Schwartzll 19891 ; lAschwandenl 2002 for reviews). A 
small fraction of the energy lost by these particles (and also 
by high-energy electrons) is emitted as gamma-ray contin- 
uum and gamma-ray lines. 

Similarly to the case of hard X-ray emission, it is dif- 
ficult to predict the gamma-ray luminosity of a cluster of 
spun-up stars. In solar flares, comparable amounts of energy 
are emitted in the hard X-ray and gamma-ray bands, which 
reflects the more fundamental fact that similar amounts of 
energy are injected into acceleration of electr ons and ions 
(see §7 and in particular Figs. 86 and 94 in lAschwandenl 
2002). Therefore, following our suggestion that a GC stellar 
cusp may produce up to ~ 10 33 erg s" 1 in hard X rays, it is 
possible that a similar luminosity is also emitted in gamma- 
rays, mainly at 0.5-10 MeV. 

This predicted flux of gamma quanta, 
<;10 _7 (MeV/S) s" 1 cm" 2 (where E is the photon en- 
ergy), is too weak to be detected by present-day gamma-ray 
telescopes. For example, the SPI spectrometer aboard the 
INTEGRAL observatory, after observing the Galactic bulge 
region for more than 10 Ms, detects 511 keV positron- 
annihilation emission with a flux ~ 10 -3 s _1 cm -2 and 
provides upper limits ~ 10 -4 s _1 cm -2 on gamma-ray emis- 



sion in the 0.511-10 MeV energy range |Churazov et al.l 
l201ll) . Moreover, the angular resolution of the existing 
gamma-ray instruments is by far insufficient to study the 
Sgr A* region. Therefore, detection of gamma-ray emission 
from a GC stellar cusp may be a matter of distant future. 

As concerns the scenario of a hot acc r etion f low onto the 
SMBH, iMahadevan. Naravan. fc Krolikl l|l997l) attempted 
to predict the gamma-ray luminosity of Sgr A* in the ADAF 
context. Some of their models resulted in a strong flux 
~ 10~ 7 s" 1 cm" 2 of ~ 100 MeV photons due to the de- 
cay of neutral pions, with a much weaker emission being 
produced at both lower and higher energies. However, we 
feel that these estimates should be regarded as indicative at 
best due to the many assumptions involved. 

4.2.2.5 Radio emission Related to the production of 
hard X rays and gamma-rays is the question of radio emis- 
sion. The radio emission of stellar coronae has been well 
studied, so we can make relatively robust estim ates in this 
case. It has been shown (|Guedel fc Bend 1 19931 ) that across 
several decades in luminosity there is a linear correlation be- 
tween the soft X-ray and radio (6 or 3.6 cm) luminosities, 
Lx/Lr ~ 10 5 ' 5 ,the same for different classes of stars. More- 
over, roughly the same relation seem s to hold for large coro- 
nal flares. Indeed. lOsten et all l|2007l ) discuss that during gi- 
ant flares of RS CVn binaries such as II Peg, V711 Tau and 
UX Ari, the radio lu minosity exceeded the quiescent levels 
by a factor of > 100 (|Waldram et al.ll2003l ; iRichards et all 
2003), i.e. by nearly the same factor as the X-ray luminosity. 

Given that the radio spectra of coronally active stars are 
approximately flat in the centimeter-d ecimeter range (e.g. 
iGarcfa-Sanchez. Paredes. fc Ribdl 120031 ) , and using the X- 
ray-radio correlation mentioned above, we may estimate the 
radio luminosity of a GC stellar cusp at z/Z/„<10 28 erg s — 1 . 
For comparison, the measured quiescent luminosity of 
Sgr A* is much higher, namely ~ 4 x 10 32 , ~ 5x 10 31 
and ~ 6 x 10 30 erg s" 1 at 3.6, 26 and 90 cm, respectively 
|An et al.ll2005l) . 

We conclude that a GC stellar cusp is unlikely to con- 
tribute significantly to the radio emission of Sgr A* and 
thus to affect the generally accepted view that this emission 
is produced in the inner accretion flow onto the SMBH. 

We summarise our comparison of the stellar cusp and 
hot accretion flow models in Table O where we indicate 
which of the observable properties of Sgr A* can possi- 
bly be accounted for by either model. We see that some 
properties (namely the iron 6.4 keV line) favour the stellar 
model, others (radio emission and strong flares of Sgr A*) 
are only consistent with the hot accretion flow model, and 
some (hard X-ray emission) allow both possibilities at the 
moment. This suggests that coronal radiation from spun-up 
stars may dominate in the extended X-ray emission around 
Sgr A*, whereas the hot accretion flow onto the SMBH possi- 
bly contributes significantly to the harder X-ray component 
and is also resposible for the phemonena spatially associated 
with the close vicinity of Sgr A*, such as the radiosource it- 
self and strong flaring activity. 

4.3 Implications for gas accretion onto the SMBH 

If a major fraction of the extended X-ray emission from 
Sgr A* is produced by a high-density stellar cluster, it will 
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Table 3. Consistency of stellar cusp and hot accretion flow mod- 
els with available data. 



Component 



Stellar cusp Accretion flow 



Quiescent emission 



X-ray radial profile 


yes 


yes 


X-ray luminosity 


yes 


yes 


Soft X-ray spectrum 


yes 


yes 


Hard X-ray spectrum 


? 


? 


Fe 6.4 keV line 


yes 


no 


Radio source 


no 


yes 


Faraday rotation limit 


yes 


? 


Variable 


emission 




Strong flares 


no 


yes 


Weak flares 


yes 


yes 



have a number of implications for the growth of the central 
SMBH, as considered below. 



4-3.1 Accretion rate 

Previously, the accretion rate through the outer boundary of 
the Bondi-type accretion flow o nto the SMBH was estimated 
at M out ~ 1O _5 M0 yr -1 (e.g. lYuan. Quataert. fe Naravanl 
2003), based on the density and temperature of the presum- 
ably diffuse gas at Rb ~ 0.05 pc, as measured with Chan- 
dra. As we already discussed in fjl] this accretion rate in 
combination with the very low (<^10 36 erg s _1 ) bolometric 
luminosity of Sgr A* puts a very strong upper limit on the 
radiative efficiency of accretion onto the SMBH and hence 
strong requirements for ADAF-type models for Sgr A*. On 
the other hand, Faraday rotation measurements in the direc- 
tion of Sgr A* imply th at accretion onto t he SMBH occurs 
at M in ^lO" 7 M yr _1 (|Bower et al.ll2003l ). 

To satisfy these observational constraints, so-called ra- 
diatively inefficient accretion flow models were suggested 
for Sgr A*. In these models, in contrast to the orig- 
inal ADAF concept, very little mass available at large 
radii actually accretes onto the black hole, while most of 
the gas outflows or circulates in convective motions (see 
lYuan. Quataert. fe Naravanl I200H and references therein). 
Hence, in this scenario, the SMBH accretes gas at a rate 
M in <;10~ 7 A'/Q yr -1 < M ou t, which allows one to both ob- 
tain the necessary X-ray luminosity and to satisfy the rota- 
tion measure limit. 

In our scenario, if the bulk of the quiescent X-ray emis- 
sion from Sgr A* is produced by spun-up stars, then the ex- 
ternal accretion rate, M out , itself may be significantly lower 
than believed before. In addition, since the stars are now re- 
sponsible for the bulk of the X-ray luminosity of Sgr A*, the 
internal accretion rate, Mi n , will be lower than thought be- 
fore. This should be taken into account and can actually be 
helpful in future developments of hot accretion flow models 
for Sgr A*. 



4-3.2 Gas supply by stars 

On the other hand, the spun-up stars themselves can pro- 
vide gas for accretion onto the SMBH via stellar winds. Al- 
though mass loss from solar-like stars increases with increas- 
ing coronal activity, this trend reverses near the saturation 



limit of coronal activity, so that the maximal loss rate of 
a few 1O~ 12 M0 yr -1 is experienced by ~ 1 Gyr old stars (in 
the solar neighbou rhood), i.e. by moderately fast rotators 
l|Wood et al.ll2005h . Therefore, since the total mass of stars 
within 0.05 pc of the SMBH is at most a few 1O 4 M , the 
integrated mass loss rate by stellar winds from the late-type 
MS stars cannot exceed 1O _7 M0 yr -1 . However, as follows 
from our preceeding discussion is £14.3.11 even with this low 
injection rate stellar winds may contribute significantly to 
the total accretion rate of the SMBH. 

A potentially more interesting mechanism of gas sup- 
ply is tidal interactions of stars. Indeed, SPH simulations 
suggest that collision-like encounters, namely those with 
R p ~ R*, not only cause a strong spin- up but also a sig- 
nificant los^_cfjnas^_i^m_the star(s), ~ O.IMq per collision 
(AK01. iFreitag fe Benzl 120051 '). For the mass density profile 
(equation [TJ with A ~ 3 x 1O 6 M0 pc -3 , such close encoun- 
ters will occur every ~ 3 x 10 4 (~ 2 x 10 4 ) years within 
0.02 (0.1) pc of Sgr A* (most collisions take place in the in- 
nermost 0.02 pc). Since these encounters occur well within 
the SMBH sphere of influence, we can expect that every 
few 10 4 years <0.1Mq of stellar debris will accrete onto the 
SMBH and produce an outburst. Less frequently, even larger 
amounts of gas (<1M ) can be provided by catastrophic 
collisions of stars, i.e. those leading to their complete dis- 
truction, and by close collisions of stars with neutron stars 
and stellar-mass black holes, which are likely to be present 
in significant numbers in the central stellar cluster. 

The injection of ~ O.l-IM© of gas from a stellar colli- 
sion onto the SMBH can cause an outburst at a significant 
fraction of its Eddington lumonisity. Indeed, the free-fall 
time to the SMBH from a typical stellar collision distance of 
~ 0.02 pc is ~ 20 years. Since the infalling gas will need to 
get rid of its angular momentum on its way to the black hole, 
the accretion episode will probably last ~ 10 2 years. This 
implies an accretion rate of ~ 1O _3 -1O _2 M0 yr -1 (if all of 
the gas actually reaches the black hole), or ~ 10~ 2 -10 _1 
of the critical accretion rate for the Ri4x 10 6 M Q SMBH. 
Therefore, stellar collisions can cause outbursts of Sgr A* 
with luminosity <10 42 ~ 43 erg s _1 every ~ 10 4 -10 5 years. 
Even stronger but shorter outbursts are expected to occur 
at a similar rate due to tidal distruction of stars in the loss 
cone by the SMBH fsee lAlexanderll2005l for review). 

The detection of reflected X-ray emission from the Sgr 
B2 molecular cloud in the GC region implies that Sgr A* ex- 
perienced an outburst with a luminosity of few x 10 39 erg s _1 
(2-200 keV) that l asted at least ~ 10 years and ended some 
100-300 years ago (jSunvacv, Markcvitc h. fe Pavlinskvlll993l ; 
iRevnivtsev et al.l 12004 iTerrier et al.ll2010r ). Such an event 
could possibly be caused by accretion onto the SMBH of gas 
torn off from a star that experienced a close collision with 
another star or stellar remnant in the nuclear cluster. More- 
over, the data also allow for the possibility that this outburst 
lasted longer than 10 years and was stronger at its outset. 
In such a case, it might have been c aused by a catastrophic 
collisi on of stars. We also note that ICramphorn fe Sunvaevl 
l|2002D searched for traces of X-ray outbursts of Sgr A* in the 
past, looking for delayed X-ray emission reflected from the 
neutral atomic and molecular gas distributed over the Milky 
Way. They concluded that there was no prolonged X-ray ac- 
tivity of Sgr A* at or above one per cent of its Eddington 
level (i.e. ~ 5 x 10 42 erg s" 1 ) in the last ~ 80,000 years. 
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This limit is consistent with the above estimates for tidal 
interaction and destruction of stars. 



5 CONCLUSION 

We have demonstrated that a significant or even dominant 
fraction of the extended X-ray emission from Sgr A* can 
be produced by a putative cusp of low-mass MS stars. To 
account for the bulk of the extended emission, there must 
be ~ (2-3) xlO 4 Sun-like stars within a column of 1" radius 
centered on Sgr A*, of which roughly every ten's star is pre- 
dicted to be rapidly rotating and producing strong thermal 
X-ray emission with kT< a few keV as a result of tidal spin- 
ups caused by close encounters with other stars and stellar 
remnants. Turning the argument around, the Chandra data 
place an interesting upper limit on the space density of (cur- 
rently unobservable) low-mass MS stars near Sgr A*. This 
bound is close to and consistent with current constraints on 
the central stellar cusp provided by NIR observations. 

The model is based on well-understood physical mech- 
anisms, namely tidal spin-up due to close fly-bys of stars 
and stellar remnants, generation of coronal activity by stel- 
lar rotation and braking of stellar rotation by magnetised 
wind. 

We have discussed potentially observable signatures of 
the stellar cusp model in comparison with the scenario of a 
hot accretion flow onto the SMBH. We demonstrated that 
the cumulative luminosity of spun-up stars in the central 
cusp can be variable by tens of per cent or even by a factor 
of few on hourly and daily timescales as a result of giant 
flares occuring on various stars. This variability could be 
searched for in archival Chandra and XMM-Newton data. 
Moreover, it is possible that some of the weaker (Lhx ~ a 
few xlO 33 erg s" 1 ), ~ 1 hour long X-ray flares from Sgr A* 
were produced by stellar coronae. There is also a chance to 
detect variable coronal emission from the GC stellar cusp at 
energies above 10 keV with the upcoming NuSTAR mission. 

We have tentatively detected a 6.4 keV line with an 
equivalent width of HOI4Q e V in the Chandra spectrum of 
Sgr A*. We argued that such a line can be expected in the 
stellar cusp model as a result of irradiation of the stellar 
photospheres by the hard X-ray coronal radiation. On the 
other hand, such emission is not expected to be produced in 
a hot accretion flow onto the SMBH. It is thus very impor- 
tant to verify the presence of a 6.4 keV line in the Sgr A* 
spectrum with more data. The planned ultradeep Chandra 
HETG observations of Sgr A* (PI: Baganoff) could solve 
this task. 

On the other hand, the stellar cusp model cannot ac- 
count for the radio emission and strong multiwavelength 
flares of Sgr A*. These are likely associated with the inner 
regions of a hot accretion flow onto the SMBH. Moreover, it 
is not yet clear if an ensemble of spun-up stars can explain 
the hard X-ray emission (at energies above ~ 4 keV) as 
easily as the softer emission from Sgr A*, since the spectra 
of typical coronally active stars may be insufficiently hard. 
We have suggested that giant stellar flares might provide 
a significant contribution to the total X-ray flux and thus 
significantly harden the cumulative spectrum of the stellar 
cusp, but this issue needs further study. 

Finally, it is plausible that spun-up stars make up the 



bulk of the softer X-ray emission from Sgr A* while a hot 
accretion flow contributes significantly at higher energies. 
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